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Abstract 

Ionization cooling in a straight beamline reduces the transverse emittance of a beam, and has Uttle effect on the 
longitudinal emittance (generally, in fact, it increases the longitudinal emittance). Once the beamline bends, the 
introduction of dispersion creates a coupling between the transverse and longitudinal planes. If this coupling is 
handled properly, one can achieve cooling in all three phase space planes. This is usually done by placing a wedge- 
shaped absorber in a region where there is dispersion. I will demonstrate using an eigenvalue analysis that there are 
other configurations of dispersion and absorber shape that will achieve ionization cooling in all phase space planes. 
In particular, I will show that a one can even achieve cooling in all phase planes with a parallel-faced absorber 
in a dispersion- free region. I will use perturbation theory to approximate the change in the cooling rates due to 
longitudinal-transverse coupling. I will then describe how the cooling of longitudinal oscillations can be understood 
via the projection of the "longitudinal" eigenmodes into the transverse plane. 

Key words: ionization cooling, coupling, emittance exchange, synchro-betatron resonance 
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1. Introduction 

Ionization cooling is a method for the rapid re- 
duction of the emittance of a beam by passing the 
beam through material (hereafter called an "ab- 
sorber") [1,2,3,4]. Its primary application has been 
to muon beams for a neutrino factory or muon col- 
lider, but it has been contemplated for other appli- 
cations as well [5,6,7]. 

In a straight beamline, ionization cooling reduces 
only the transverse emittance of a beam, generally 
having little effect on the longitudinal emittance 
(in fact, generally making it somewhat worse). It 
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reduces the transverse emittance because the to- 
tal momentum, including the transverse, is reduced 
when when the beam passes through material, but 
when the momentum is restored by an RF cavity, the 
transverse momentum is left unchanged, and thus 
reduced from its value before the absorber. It is de- 
sirable, especially for a muon collider, to reduce the 
longitudinal emittance as well. 

There is a well-known method for accomplish- 
ing this, often referred to as "emittance ex- 
change" [4,8,9]. This is generally accomplished by 
using a wedge-shaped absorber in a region with 
non-zero dispersion, as shown in Fig. 1. Particles 
with higher energy will pass through more material 
than particles with lower energy as a result of the 
dispersion. Thus, the energy spread in the beam 
will be reduced. There is a cost to this, however. 
Imagine that all the energy spread were removed by 
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Fig. 1. Dispersion generated by a bending magnet (left) 
means that particles with different energies have different 
positions. The wedge-shaped absorber (right) then reduces 
the energies of higher energy particles more than it reduces 
the energies of lower energy particles. 

this. Then there would be no way of removing the 
spread in transverse position that was introduced 
by the dispersion. Thus, reduction in longitudinal 
omittance comes at a cost: the horizontal emittance 
is increased. If done properly, the cooling accom- 
plished transversely by an absorber can be shared 
between the longitudinal and transverse planes, giv- 
ing emittance reduction in all phase space planes. 

This paper will demonstrate that achieving cool- 
ing in the longitudinal plane can be accomplished by 
a more general class of methods than just having dis- 
persion in transverse position at a wedge-shaped ab- 
sorber. It will demonstrate that cooling in all planes 
happens when there is an appropriate coupling be- 
tween the longitudinal and transverse planes. The 
paper will explore multiple ways of generating that 
coupling and compare their effectiveness. 

I will first describe a simplified lattice that will 
be used in this paper. I will then ascertain whether 
the lattice is cooling in all phase space planes by 
examining the eigenvalues for the linear matrix rep- 
resenting that lattice. I will compute the eigenval- 
ues approximately using perturbation theory, where 
the small parameter in the expansion is the degree 
of coupling between the longitudinal and transverse 
planes. I will study coupling that is generated by dis- 
persion and an appropriately shaped absorber, and 
coupling generated by dispersion in the RF cavity. 
I will then look at the eigenvalues for cases when 
the perturbation expansion might not be accurate. 
Finally, I will develop a more physically intuitive 
understanding of how transverse cooling can reduce 
the longitudinal emittance by examining the motion 
in the coupled system using the eigenvectors of the 
matrix that represents the lattice, which will give a 
picture of the coupled motion when it is projected 
into the longitudinal and transverse planes. 

This paper does not address other issues of cool- 
ing lattice design, in particular the final equilibrium 
emittance which results because of multiple scat- 



tering and energy straggling. The paper is intended 
only to present a wider range of ideas on how cool- 
ing can be accomplished in all phase space planes, 
not to present a full generalized theory of ionization 
cooling or a working cooling lattice design. 



2. Lattice 

The cooling lattice will consist of a sequence of 
identical cells with four sections: 

(i) An absorber, which reduces the energy of the 
particles 

(ii) An RF cavity, which restores the energy lost 
in the absorber and provides longitudinal fo- 
cusing 

(iii) A group of magnets transporting the beam 
from the absorber to the next RF cavity 

(iv) A group of magnets transporting the beam 
from the cavity to the next absorber 

The details of what is in the "group of magnets" will 
not be important for this discussion. 

There is a planar reference curve which defines 
the coordinate system that we are using. The ver- 
tical coordinate y is perpendicular to the plane in 
which the curve lies. The horizontal coordinate x is 
perpendicular to the vertical and to the tangent to 
the curve, and stays on the same side of the curve. 
The arc length along the curve, s, is the indepen- 
dent variable for the equations of motion. For the 
examples in this paper, the vertical magnetic field 
at the reference curve Bo{s) will be such that a par- 
ticle with a reference energy £"1(5) which starts out 
on the curve and moving tangent to the curve con- 
tinues to do so. This particle will be referred to as 
the reference particle. Note that the reference en- 
ergy depends on the position in the lattice, due to 
the absorber and the RF cavity. 

The magnetic field will be midplane symmetric, 
meaning that the vertical field is symmetric under 
y — !■ —y, and the other magnetic fields are antisym- 
metric under y —y. There are two consequences 
of this. The first is that this is not a solenoid focused 
lattice, which is typical for an ionization cooling lat- 
tice, but a quadrupole-focused lattice. The second is 
that the vertical motion is, to lowest order, decou- 
pled from the horizontal and longitudinal motions. 
This second reason is the fundamental purpose in 
choosing a midplane symmetric lattice. The quali- 
tative results found here will continue to hold for a 
solenoid-focused lattice. 
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The reference particle will lose an energy AE in 
the absorber and will gain back that same amount 
energy in the RF cavity. The energies in the mag- 
nets between these elements will be E± = Eq ± 
AE/2, the corresponding total momenta are p± = 

\J E'j^/c^ — {racY, and po = ^/^a/'^^ ^ {mc)^. In 
these expressions m is the particle mass and c is the 
speed of light. 

2.1. Matrices for Lattice Sections 

I want to determine whether for small deviations 
from the particle with coordinates zero, the ampli- 
tude of oscillations is growing or reducing. Thus, this 
paper will only keep results to first order in these 
deviations, and will use matrices to represent how 
these deviations change while propagating through 
the lattice. The lattice is a sequence of identical cells, 
and the eigenvalues for the matrix for an entire cell 
will determine the characteristics of the motion. 

I will denote the absorber location by a subscript 
"a," and the cavity location by a subscript "c." The 
matrix describing the motion through the magnets 
from the absorber to the cavity is Aica, and the 
matrix from the cavity to the absorber is A^ac- 1 will 
assume a kind of reflection symmetry at the absorber 
and cavity such that these matrices can be written 
as 



A^ac 



(1) 
(2) 
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These quantities are related to the usual accel- 
erator quantities by bi = j3x{si) /p^{si)^ df = 
D{si)/[pr{s^)pr{si)c] and df = D'{s ,)/[PM)c], 
where pr{s) = ^E^{s)/c'^ - {mcY, pr{s) = 
Pr{s)c/ E^{s), 0x{s) is the Courant-Snydcr beta 
function, and D{s) is the dispersion. Note that 
since Pr(sa) = P± depending on which side of the 
absorber one is on, the beta fimctions on each side 
of the absorber are actually slightly different. The 
same type of difference exists for the cavity and the 
dispersion. 

Assuming that the only effect of the RF cavity 
and absorber is to shift the energy of all particles by 
ztAE, the matrix for the full cell from the absorber 
back to the next absorber is 



V^BJi{n-+n+,R 



56 



(6) 



Thus, the transverse phase advance per cell is /i_ -I- 
and R^Q -\- R^q is approximately r]cTE/{pc)^, 
where rjc is the frequency slip factor, and T is the 
time for the reference particle to go through one cell. 
The momentum p and energy E are intentionally 
left ambiguous: rjc is only defined for a fixed energy. 
The RF cavity is described by a matrix 



V 



10 
1 
1 

coVcosd) 1 



(7) 



where V is the maximum energy gain in the cavity, 
the RF phase is (f>, where (/) = is the phase for zero 
energy gain, and co is 27r times the RF frequency. 

The absorber reduces the energy of the particles, 
maintaining their direction. The horizontal position 
will thus be the same as it would be if the absorber 
were a drift. The time to traverse the absorber is 
well approximated by assuming that the momentum 
is p+ to the center of the absorber, then p_ there- 
after. Thus, it is a good approximation to treat the 
absorber as a thin clement, with surrounding drifts 
absorbed into Mac and A^ca- 

Since the particle direction is maintained, the 
transverse momentum will be reduced in the ab- 
sorber by the same factor that the total momentum 
is reduced, = p-/p+. Particles whose energy dif- 
fers from the reference energy will receive a slightly 
different energy loss than the reference particle, 
since the energy loss per unit length in the absorber, 
dE/dx, depends on energy. To linear order, the en- 
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Fig. 2. Geometry of the absorber, sincji = px/v+^ where Px 
is the horizontal momentum before the absorber. 

ergy loss of a particle with energy E just before the 
absorber is 



f-1 




(dE\ 


I dx J 




\ dx I 



= AE + {1 - K,){E - E+) (8) 

Furthermore, if the faces of the absorber are not 
perpendicular to the reference orbit, there may be an 
energy loss which depends on the horizontal position 
or angle (again, I'm ignoring the vertical motion). 
Using the geometry from Fig. 2, the change in the 
energy is a+a; + a-Lpx/{2p+) [10], where 

dE 

a± = — — (tan 6*2 ± tan 6*1 ) . (9) 

dx 

These define the matrix for the absorber, which I 
will denote A: 



A = 



1 




























1 






a+ a_L/(2p+) Kz 



(10) 



3. Eigenvalue Analysis 

The transfer matrix for the full lattice,including 
the absorbers and the RF cavities, is 



I will be analyzing the eigenvalues of this matrix to 
determine if I am achieving cooling in both planes. 
First of all, it is clear that the product of the eigen- 
values and that the eigenvalues will come in 
complex conjugate pairs. 

Begin with a basic lattice, with no dispersion, and 
with the absorber face angle 9i = 62 = 0- The char- 
acteristic polynomial for the matrix for the cell is 

(A^ - A(l + + {R^,. + K,R+.)iuVcos(j)) + Kz) 

(A^ - A(l -|-Kx)cos(/i+ +KX)- (12) 

The roots of the characteristic polynomial are the 
eigenvalues. 

Prom this polynomial, we can see that the eigen- 
values will be of the form 



K,e 



(13) 



It turns out that the eigenvalues remain the same 
if d^ and d^ are nonzero, as long as one minor cor- 
rection is applied: in A4ac, d^ should be replaced 
with dav^' -M-ca., d^ should be replaced with 

Thus, ^/K^ is the magnitude of the transverse 

eigenvalue, and ^/k^ is the magnitude of the longitu- 
dinal eigenvalue. Clearly < 1; unfortunately, in 
most cases of interest for cooling, > 1. However, 
the product Ka;^^ < 1. > 1 because for low en- 
ergies, {d/dE){dE/dx) < 0. While {d/dE){dE/dx) 
becomes positive for higher energies, the relative en- 
ergy loss for a given absorber is less at higher ener- 
gies, making Kx closer to 1, and the product Kx^z 
closer to 1. 

Since the product k^k^ < 1, if one could gener- 
ate coupling between the horizontal and vertical mo- 
tion, it would be possible to reduce the magnitude 
of the "longitudinal" eigenvalues (they begin to lose 
their identity with coupling), at the cost of raising 
the magnitude of the "transverse" eigenvalues. One 
could do so in a way that made the magnitudes of 
both sets of eigenvalues less than 1, and thus gave 
cooling in all phase space planes. 

If the characteristic polynomial for a matrix is 
Eq. (12) plus some additional polynomial /(A), then 
one can compute the change in the eigenvalues to 
lowest order. The lowest-order change in the magni- 
tude of \x = J^x 
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and the lowest-order change in the magnitude of 
K = ^J^Tze^^' is 



4 



CSC /i2 



Im 



(A. - A,)(A, - A*) 



(15) 



Since a change in the magnitude of \z leads to a 
corresponding change in the magnitude of A^, one 
only need examine the change in one eigenvalue. 

Note that when A^ is close to A^ or A*, there will be 
a large change in the eigenvalues. This corresponds 
to a linear coupling resonance between the longitu- 
dinal and horizontal. For our perturbation analysis, 
we will assume that we are sufHciently far from that 
resonance condition, but this will be of interest later 
in the paper nonetheless. 

Since / is the difference between the characteris- 
tic polynomial for the full system and Eq. (12). wo 
know something more about its properties. First of 
all, it is a third order polynomial, since all character- 
istic polynomials of a 4 x 4 matrix are fourth order 
polynomials with leading order term A^. Secondly, 
the constant term in the characteristic polynomial is 
the determinant, and the determinant of the matrix 
is KxKz in all cases. Thus, / will have no constant 
term. We can thus write /(A) as /aA^ + /2A^ + /lA. 

I will examine the change in the magnitude of X^. 
Using the polynomial expression for /(A), it is 



"i.^UxKziUzh - /l) COS IJLx 
- 2Kz{Kxh - h) COS Hz + \/l^{Hz - Hx)h 
'^[k.z + K.X- 2V«a:Kz COs(/i^ - /X^)] 

[Kz + Ka: - COs(/Ua; + Hz)] 



(16) 



There arc two methods for generating coupling 
between the phase space planes. The most common 
method is to create an absorber with nonzero values 
for the angles 6i and 62 (in most cases with 61 ^02, 
thus referred to as a "wedge"). Another method is 
to generate the coupling by having dispersion in the 
RF cavity, which docs not require any change in the 
shape of the absorber. 



3.1. Rotated Absorber Faces 

The additional terms in the characteristic poly- 
nomial when ^0 or Oi are nonzero but the dispersion 
at the cavity is zero are 



—=[X'{COSIIxO - KxS+) 
V l^x 

+ A^ cosiixo{kIs+ - S-) 

+ XKx{S- — Kx COS/Xa:o)] 

-I- a+d^ba sin HxoiX^ - A^(kxS+ -|- S-) -h A/t^] 



2p_^_y/K^ba 



sin/Xj,o[A - A {kxS+ + s_) + Xk^] 



+ -T, — ^ A'^(cosM^o - s+) 
2p+ 

+ X'^ cosiJxo{s+ - S-) 

+ X{s- - coshm)], (17) 

where ^xo = M+ + ^^'^ s± = 1 + R^qUjV cos cj). 
Note that from Eq. (12), KzS+ + s_ = 2yKjcos /i^. 
The numerator of Eq. (16) for this case is times 

^^^[s_ + KxS+ - {1 + rtx) COS/Xxo] 
V f^x 

[k^{S- + S+K^) - {kI + K^) COS/ixo] 

- a+(i^6a(K2 - Kx) sin fiM 

[.S_ + Ka..S+ - {1 + Kx) COS/Xa;o] 

- Kx ) sm Hxo 



S- + KxS+ - {1 + Kx) COS/ia:o] 

a-LdP, , , 

+ -5 [S- + KxS+ -{1 + Kx) COSfJ,x0\ 

[S-+S+K;,-{l + Kz)cOSIJ.x0]- (18) 



Since « 1 in most cases, the quantity 

be treated as small. Thus, the second 
and third terms in Eq. (18) can be neglected, as 
long as one is far from the linear synchro-betatron 
resonance. Equation (18) indicates that if one has 
a dispersion in position but not momentum at the 
absorber, the absorber should have 61 = 62 (i.e., 
a wedge shape), and d\ positive if the dispersion 
is positive. This is the usual method of "emittance 
exchange." 

If, on the other hand, there is no dispersion in po- 
sition but there is dispersion in horizontal momen- 
tum, one slioidd instead use a rotated slab (i.e., 0\ = 
—^2)- Particles with larger momentum will have an 
angle in one direction; if the absorber is rotated so 
there is more material along that direction, there 
will be a greater energy loss for those particles. 

One might ask whether it is easier to construct 
a lattice that makes the first term large than it is 
to construct a lattice that makes the fourth term 
large. In fact, either term can be made large, and 
which one is used depends on the desired properties 
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of the lattice. Take, for instance, the RFOFO cool- 
ing ring described in [11]. {dE / dx)d^ for that lat- 
tice is about 0.027. If one finds the largest value of 
the momentum dispersion, then {dE/dx)Ld^/ (2p+) 
is around 0.002. While it may seem that the for- 
mer is significantly larger than the latter, the fact 
is that the lattice was designed with a ring shape, 
which tends to generate position dispersion rather 
than momentum dispersion. One can construct lat- 
tices where {dE / dx)Ld^/ {2p+) is larger, especially 
when the beta function at the absorber is small, by 
having one half of the lattice cell bending in one di- 
rection and and the other half bending in the reverse 
direction. 

3.2. Dispersion in the RF Cavity 
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Fig. 3. Magnitude of matrix eigenvalues as a function of 
(//+ + n- ) / (27r) for a lattice with dispersion at the RF cav- 
ities and an absorber with 0i = 62 = 0. 



If instead, 9o = Oi =0, but we have dispersion at 
the cavity, there is still coupling between the longi- 
tudinal and transverse planes, in this case generated 
by the dispersion at the RF cavity. The additional 
terms in the characteristic polynomial are 

— — A(A- 1)(A- K^) 

{[W)' + fe?ra'](«x + l)sin(M++M-) 
-[(d-)2-62(dP)2](«,-l)sin(M+-M_) 

- 26cd?(iP(K^ - 1) cos(/x+ - II.)} (19) 

The numerator of Eq. (16) is then 

ojV cos ' 



{Kz - Ki;)[2k2 COSyU^ - V«r(l + K^) 

blidP,f]{n, + l)sm{f,+ +f,_) 
■^c(rf?)'](«.-l)sin(M+-M-) 



2bc 

- Mr 

- 2hcdldl{K^ - 1) cos(Ai+ - ii-)} 
For the RFOFO ring in [11], the quantity 
cos 



26c 



(20) 



(21) 



is around 1.7 x 10~^. The product of the next two 
factors in Eq. (20) is around —0.02. Thus, compared 
to an absorber with substantial rotations in its faces, 
the strength of the coupling generated by dispersion 
in the RF cavities is relatively weak. The RFOFO 
ring was not necessarily designed for maximizing 
coupling in this way, but it is unlikely to be possible 
to increase the dispersion sufficiently while keeping 
a reasonably small beta function in the lattice. One 
may be able to increase the RF voltage somewhat, 
but certainly not enough to change the overall re- 
sult. 



3.3. Running Closer to Resonance 

The perturbation theory analysis assumed that 
one was far from the point where any two of the 
eigenvalues of the matrix are close to each other. If 
they do become close, then the estimates above are 
not valid, and the eigenvalues should be computed 
directly. 

Begin with the case where there is dispersion in 
the RF cavity and = ^1 =0. Figure 3 shows the 
magnitude of the eigenvalues as a function of /i+ + 
where I am taking /i+ = Note that when 
(/j,+ + /i_)/(27r) « 0.08, which is the synchrotron 
tune, then the magnitudes of all the eigenvalues be- 
come less than 1 [12]. What one is seeing here is 
the manifestation of a c;oupling resonance in a non- 
conservative dynamical system. Note that the same 
resonance phenomenon at (/i+-|-/z_)/(27r) « 1—0.08 
results in the longitudinal eigenvalue becoming more 
unstable: instead of the resonance pulling the eigen- 
values together, it pushes them further apart. The 
signs of these effects are consistent with what one ex- 
pects from the perturbation calculation in Eq. (20): 
the sign of sin(/i_|_ + H-) gives the dominant effect 
for determining whether the magnitude of the longi- 
tudinal eigenvalue increases or increases. The loops 
in Fig. 3 near the tunes of and 0.5 are not impor- 
tant for this discussion: they are simply the linear 
resonances of the transverse lattice. 

This example is furthermore an example of how a 
difference resonance {vx — Vz = k'm this case) tends 
to lead to often innocuous coupling between planes, 
whereas a sum resonance {v^ + i^z = k in this case) 
can lead to unbounded growth. The coupling be- 
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Fig. 4. Magnitude of matrix eigenvalues as a function of 
+ ;u_)/(27r) for a lattice with dispersion at a wedge 
shaped absorber. 

tween the planes is in fact an advantage in this case, 
since one is trying to make the transverse damping 
affect the longitudinal plane. 

Having a low horizontal tune is generally imprac- 
tical for cooling channels, since that tends to give a 
larger beta function at the absorber, which results 
in a larger equilibrium emittance due to multiple 
scattering [3,4]. However, one could instead run in 
the passband with horizontal tunes from 1 to 1.5, 
which would both create a low beta function and 
ensure that coupling pulled the magnitudes of the 
eigenvalues closer to each other. It would also allow 
for a longer cell, making it possible to increase the 
synchrotron tune, putting it further from the inte- 
ger resonance. It may be more difficult to have a 
large momentum acceptance in such a lattice, how- 
ever, since the passband from tunes of 1 to 1.5 gen- 
erally would have a smaller relative momentum ac- 
ceptance than the lower passbands. 

This example is meant more as a proof of prin- 
ciple, demonstrating that one can in principle use 
any method of longitudinal-transverse coupling 
to achieve cooling in all degrees of freedom. The 
method may, however, have practical applications 
in cases where it is impractical to control the shape 
of the absorber. One example might be the use of a 
lithium lens in the final stages of cooling for a muon 
collider. The method may be more interesting at 
later stages of a muon collider in any case, since the 
limited width of the resonance will likely translate 
into a limited momentum acceptance for such a 
method, requiring a beam which has already had 
its longitudinal emittance reduced significantly. 

It is interesting to compare Fig. 3 with what one 
would obtain when using a wedge-shaped absorber 



with dispersion. This is shown in Fig. 4. Note that 
the eigenvalue magnitudes are closer to each other 
over the entire range of horizontal tunes (their val- 
ues without dispersion at the wedges can be de- 
duced from Fig. 3). At both the sum and differ- 
ence resonances, the magnitudes are pushed fur- 
ther apart. The difference in the behavior near the 
synchro-betatron resonances between this case and 
that shown in Fig. 3 is a result of the non-symplectic 
nature of the coupling that is generated by the wedge 
with dispersion. Eqs. (16) and (18) do an extremely 
good job of predicting the behavior in Fig. 4, in- 
cluding the resonant behavior; this is to be expected 
since the eigenvalues never get too close together, 
even near the resonance. 

4. Physical Explanation 

To get some understanding of what is going on 
physically, one should examine the eigenvectors. The 
real and imaginary parts of the eigenvectors define 
an ellipse that the particles move on. Of course, in 
the case of cooling, the radius of the particles are 
decreasing (or sometimes increasing) , but that piece 
can be taken out (i.e., the difference of the magni- 
tude of the eigenvalue from 1). To see the ellipse, 
take the real (or imaginary) part of the eigenvec- 
tor, multiply the vector by the matrix for the lattice 
cell, and divide the result by the magnitude of the 
corresponding eigenvalue. Repeat the process, and 
one will trace out an ellipse in phase space. Alterna- 
tively, if V is the eigenvector, then the ellipse is the 
set of points obtained from 

Re{t;} cos u + Im{t>} sin u (22) 

by varying u from to 27r. 

In the case where there is no coupling, the eigc;n- 
vectors have components either entirely in the hor- 
izontal plane, or entirely in the longitudinal (time- 
energy) plane. However, once one introduces any 
type of longitudinal-transverse coupling, all of the 
eigenvectors will have components in both planes. 
This is the key to what allows the reduction in trans- 
verse momentum spread that the absorber accom- 
plishes to reduce the amplitude of what would oth- 
erwise be longitudinal motion. The only effect that 
actually reduces the beam emittance is the reduc- 
tion of transverse momentum in the absorber. 

Figure 5 shows the projection of the longitudi- 
nal eigenvector (identified by the phase of its corre- 
sponding eigenvalue) at the absorber into the trans- 
verse plane. The figure shows two cases: one where 
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Fig. 5. Projection of the longitudinal eigenvector at the ab- 
sorber into the transverse plane when there is dispersion at 
the absorber. The solid ellipse (which has no area, so it ap- 
pears as a line) is for a parallel-faced absorber {0i =82 = 0), 
whereas the dashed ellipse is for a wedge-shaped absorber. 

there is dispersion at the absorber, but the absorber 
has parallel faces with 9i = 62 =0, and one where 
the absorber is wedge shaped, with 9i = 62 > 0. 
For the parallel-faced case, there is no reduction in 
the magnitude of the longitudinal eigenvalue, de- 
spite the fact that Fig. 5 shows that there is a non- 
zero projection of the longitudinal eigenvalue onto 
the horizontal plane (and horizontal momentum in 
particular). The reason that there is no reduction 
in the magnitude of the longitudinal eigenvalue is 
that the ellipse projected into the horizontal plane 
has no area. The absorber makes positive momenta 
more negative and negative momenta more positive, 
reducing the area of an ellipse. But if the ellipse has 
110 area, no reduction can occur. Thus, one needs 
the projected ellipse to have a nonzero area to see 
an effect. This is what making the absorber wedge- 
shaped accomplishes, as can be seen in Fig. 5. 

If one changes the horizontal tune to around 0.077, 
one can see from Fig. 4 that the longitudinal and 
transverse eigenvalues return almost to their values 
without coupling. This is reflected in the eigenvec- 
tors in that the area of the horizontal projection of 
the ellipse relative to the area of the longitudinal 
projection is less than it was for horizontal tunes 
away from the longitudinal tune (which was shown 
in Fig. 5). 

Figure 6 shows the ellipses for the case where 
there is dispersion at the RF cavity and not the 
absorber (this corresponds to Fig. 3). The ellipses 
are smaller when one is away from resonance, and 
larger when one is closer to resonance, reflecting the 
stronger coupling between the eigenvalues near the 



Fig. 6. Projection of the longitudinal eigenvector at the ab- 
sorber into the transverse plane when there is a parallel-faced 
absorber and dispersion only at the RF cavity. The solid el- 
lipse is for a horizontal tune of 0.75, the dotted ellipse for a 
horizontal tune of 0.09 (near one coupling resonance), and 
the dashed ellipse for a horizontal tune of 0.91 (near the 
other coupling resonance). 

resonance, in particular that the magnitude of the 
longitudinal eigenvalue is reduced (or increased for 
the higher-tune resonance) there. What the flgure 
does not illustrate is why near the low tune reso- 
nance (0.09 for the example in Fig. 6), the magnitude 
of the longitudinal eigenvalue is reduced, whereas 
near the high tune resonance, the magnitude of the 
longitudinal eigenvalue is increased. 

To understand this, one must first think about 
what area means for an ellipse in four- dimensional 
phase space. Hamilton's equations of motion can be 
written as 

^ = -JVHiz,s), (23) 

where for the example here where only horizontal 
and longitudinal dimensions arc considered, 

10 

-10 
J = (24) 

-1 
10 

One can define the area of a four- dimensional ellipse 
which is described by Eq. (22) to be tt times 

(Im{t;})^J(Re{t;}). (25) 

This gives the expected answer when the ellipse lies 
entirely in either the horizontal or longitudinal phase 
space plane, and is thus the obvious extension of 
the concept. In fact, it is really the sum of the pro- 
jected areas in the two planes. Furthermore, it is in- 
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variant under symplectic transformations of v, and 
thus reflects the area-preserving nature of symplec- 
tic transformations. 

Note that this defined, has a sign. Since 

the eigenvectors come in complex-conjugate pairs, 
there are always two areas with the same magnitude 
and opposite signs generated by the eigenvectors of 
a matrix, assuming the eigenvectors were not nor- 
malized differently. Furthermore, since the area is 
really a sum of the horizontal area and the longitu- 
dinal area, these projected areas individually have 
signs as well. It turns out that when the longitudi- 
nal eigenvalue gets closer in magnitude to the trans- 
verse one, the areas have the same sign, whereas 
when they are pushed away from each other, the ar- 
eas have opposite signs. 

Physically, this is because the sign of the ellipse 
area is related to the direction which a particle moves 
around the ellipse in phase space, as can be seen from 
Eq. (22). There is a "physical" direction that parti- 
cles move around an ellipse, which is reflected in the 
signs of the nonzero elements in J. Normally, par- 
ticles with positive horizontal momentum will have 
an increasing horizontal position; particles with a 
larger energy will have a shorter time of flight (ig- 
noring momentum compaction). When the signs of 
the projected areas (horizontal and longitudinal) are 
identical, the particles are moving in the "physical" 
direction on both the horizontal and longitudinal 
ellipses. If, on the other hand, the projected areas 
have opposite signs, then longitudinal motion in the 
"physical" direction will appear to have "unphysi- 
cal" motion in the horizontal plane. Thus, in par- 
ticular, the reduction in horizontal momentum that 
the absorber accomplishes will not be translated 
physically into the area reduction that one hopes 
for, since the change in horizontal position goes the 
wrong way. The result is the eigenvalue magnitudes 
getting further apart, as seen on the right hand side 
of Fig. 3. 

5. Conclusion 

I have shown that one can obtain ionization cool- 
ing in all phase space planes through two differ- 
ent methods: by generating dispersion at the ab- 
sorber and shaping the absorber appropriately (this 
is the usual method), and by using the longitudinal- 
transverse coupling resonance with a parallel-faced 
absorber. I have approximated the effect by pertur- 
bation theory, and using that have demonstrated 



that one should shape the absorber in conformity to 
the dispersion at the absorber, in particular taking 
into account whether the dispersion is in position or 
transverse momentum. I have demonstrated that it 
is the fact that the longitudinal motion has a com- 
ponent in the transverse plane that allows the ab- 
sorber to cause cooling in the longitudinal plane as 
well as the transverse. The strength of the effect is 
related to the area, including sign, of the "longitu- 
dinal" ellipse projected into the horizontal plane. 
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